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Mass spectra of ground state hadrons containing u-, d-, s-, c- quarks as well as 
some lightest hadrons containing fr-quarks are calculated on the basis of a slightly 
, modified bag model. The center-of-mass motion corrections are incorporated using a 

£> \ wavepacket projection with Gaussian parametrization of the distribution amplitude. 

We use running coupling constant and also allow the effective quark mass to be 
scale-dependent. The impact of these modifications on the hadron mass spectrum 
is investigated. A comparison of the predicted mass values with the experimental 
data demonstrates that the modified bag model is sufficiently flexible to provide a 
satisfactory description of light and heavy hadrons (mesons and baryons) in a single 
consistent framework. 
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1 Introduction 



Over the last decade a lot of progress has been made in the experimental 
spectroscopy of heavy hadrons. Accumulation of the high statistics data by 
various experiments led to the discovery of many new states. Among others, 
even rather exotic state of two different heavy quarks (B c meson) [1,2] has been 
observed. In addition, spectroscopy of heavy hadrons serves as an important 
field to test various QCD-inspired phenomenological models of hadron struc- 
ture. One of such models is the MIT (Massachusetts Institute of Technology) 
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bag model [3,4]. There are several excellent reviews on this subject available 
[5,6,7,8,9], where one can find more information concerning basic equations, 
applications, and further developments of the bag model. 

After the first success in describing the static properties of the light hadrons 
[10], a straightforward application of the bag model to the heavy quark states 
[11] led to a surprisingly strong disagreement with the experimental data. 
Early attempts [12,13] to improve the model were of limited success. Discrep- 
ancies seemed to be of qualitative character, so one could conclude that some 
more radical modifications of the model were necessary. It was soon realized 
that the bag model was afflicted by the well-known center-of-mass motion 
(c.m.m.) problem. A part of the hadron energy calculated in the ordinary bag 
model is spurious and, consequently, the model must be corrected in some 
fashion. Such correction may lead to the substantial changes in the predicted 
mass values of the light hadrons [14,15,16]. To the best of our knowledge, 
at present there is no unambiguous method to deal with this problem. Nev- 
ertheless, approximate schemes have been widely used in various bag model 
calculations [17,18,19,20,21,22,23,24,25,26,27]. 

For the hadrons containing one heavy quark an elegant way to eliminate the 
center-of-mass motion has been proposed [28,29] (for further developments see 
[30,31]). In that approach the heavy quark occupies the center of the bag and 
the light quarks move in the colour field set up by this heavy quark. A simple 
physical picture is an attractive feature of this prescription. However, its appli- 
cability is restricted to hydrogen-like systems. Therefore, if we want to have a 
unified description of the hadrons, we need a more universal tool to deal with 
the c.m.m. problem. Although there is some controversy on this subject, we 
have chosen to follow the technique adopted in Refs. [24,25,26,27]. The essence 
of this method is to replace the bag state with the wave packet (a superposi- 
tion of plane- wave states). A similar approach to correct for the c.m.m. was 
used within the framework of the relativistic potential model [32,33,34,35]. 

The aim of this paper is to provide a unified description of the light and 
the heavy hadrons in the framework of the bag model. Besides the c.m.m. 
correction we will incorporate two other QCD-inspired improvements of the 
model: the running (i.e. scale-dependent) effective coupling constant, and the 
scale-dependent effective quark mass. The influence of these improvements on 
the hadron mass spectrum will be investigated. 

This paper is organized as follows. In the next section a modified bag model 
is described. Our results on the calculated hadron spectrum are presented in 
Section 3 along with a discussion of the influence of the modifications upon the 
bag model predictions. Some other static parameters of light hadrons (mag- 
netic moments, axial-vector coupling constant 5% and charge radii) for which 
experimental data exist are calculated and presented in Section 4. Finally, we 
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summarize our conclusions in Section 5. 



2 The model 



The bag model enables us to calculate the static properties of hadrons by 
making a number of simplifying assumptions. Usually it is assumed that the 
quarks are confined in the sphere of fixed radius R, within which they obey 
the free Dirac equation (static spherical cavity approximation). The energy of 
a hadron is given by 

4-7T 

E = —BR 3 + J2n i e l + AE. (1) 
6 i 



The first term on the right-hand side of the above equation is the bag volume 
energy that guarantees the quark confinement in the finite region, R stands 
for the radius of the bag, and B is the bag constant. The second term is the 
"kinetic" energy of quarks, rij is the number of quarks of z-th flavour, £j - the 
eigenenergy of a quark in the cavity. The last term represents the interaction 
energy of the quarks in the Abelian approximation to QCD. Minimization of 
the energy determines the bag radius R of the hadron under consideration. 

It is useful to divide AE into two parts: 

AE = E m + E e . (2) 



One, 



E m = a r 



auMu + a>ij M K 



(3) 



is the colour-magnetic part, and another, 



E e = a r 



i j>i 



(4) 



is the colour-electric (Coulomb) part of the interaction. In Eqs. (3) and (4) 
a c is the coupling constant and the sum runs over the flavour indices. For the 
benefit of the reader, below we present the expressions (3) and (4) in more 
detail, omitting tedious derivation procedures. Functions Mij(R) and Iij(R) 



3 



can be written in the form 

R 



M iJ (R) = -Jdr^ l (r)A J (r,R), (5) 



4 

o 

R 



I t ,(R) = -Jdrp' t (r)V,(r,R). (6) 
o 



Here 



m = -yPi(r)Q<(r) (7) 



is the scalar magnetization density of an i-th quark. The semiclassical vector 
potential generated by the z'-th quark has the form [36] 

Mr,R) = ^ 1 + ^§ + M l (r,R), (8) 



where 



r 



r 

J dx^(x), (9) 



o 

R 



Mi{r,R) = I dx^-. (10) 

x 6 



In Eq. (6) 

A(r) = P?{r) + QUr) (11) 
is the charge density of the z'-th quark, and 

r 

/ dxfttx), (12) 

Ft 

K(r, R) = Pl {r) (I - 1) + j dx^l. (13) 



r 

p,(r) 

o 



r 
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Pj(r) and Qi{r) in Eqs. (7) and (11) are the large and small radial functions 
of the two-component spherical spinor normalized as 



R 



Jdr [P? (r) + Q* (r)] = 1, 



(14) 



and obeying the linear boundary condition 
P t (R) = -Qi{R) 



(15) 



at the bag surface. 

In order to avoid any possible complications we have used the confined Coulomb 
Green's function [37] in the derivation of the expression (6). As a consequence, 
the value of the colour scalar potential Vi(r, R) is zero at the surface of the 
cavity: 



V i (R,R)=0. 



(16) 



The coefficients a^, fi, and f\j that specify the interaction energy of hadrons 
in Eqs. (3) and (4) can be readily calculated using the technique described in 
Ref. [10]. Parameters /; that specify the colour-electrostatic interaction energy 
between the quarks of the same flavour are 



fi = -X- m{ni - l)/2, 



(17) 



and parameters (i ^ j) are given by 

fij A Tl-iTlj, 



(18) 



where 



A 



1 for baryons, 

2 for mesons. 



(19) 



Parameters that specify the colour-magnetostatic interaction energy for 
mesons with the total spin J are 



-6 (J = 0), 
2 (J=l). 



(20) 
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Table 1 

Parameters that specify the colour-magnetic interaction energy of baryons consisting 
of the light (I = u,d), strange (s), and charmed (c) quarks. 



J 


Particle 


Quark content 


an 


a ls 


a lc 








1/2 


iV 


III 


-3 












1/2 


A 


\ / Ulbtsl 


-3 












1/2 


£ 




1 


-4 










1/2 




Iss 




-4 




1 






1/2 


A+ 

C 


L- \ lilt ) fiwhi 
\ / til vol 


-3 












1/2 


s c 


L- 1 L L to iTn 


1 




-4 








1/2 




\ / <_£ / ttfc 




-3 










1/2 


r-n/ 

1 'c 






1 


-2 




-2 




1/2 


n° r 

c 


CSS 








1 


-4 




1/2 




Ice 






-4 






1 


1/2 


cc 


sec 










-4 


1 


3/2 


A 


III 


3 












3/2 


S* 


sll 


1 


2 










3/2 




Iss 




2 




1 






3/2 


n- 


sss 








3 






3/2 


c 


ell 


1 




2 








3/2 


"* 


els 




1 


1 




1 




3/2 




CSS 








1 


2 




3/2 


"* 


lec 






2 






1 


3/2 




sec 










2 


1 


3/2 


"ccc 


ccc 
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For the baryons consisting of u-, d-, s-, and c-quarks these parameters are 
given in Table 1. For the baryons containing 6-quarks the corresponding pa- 
rameters can be easily defined by means of simple substitutions (e.g., c — > b). 
In the case of the light hadrons one can also find the parameters in Table 2 
of Ref. [10]. 

Now we describe the salient differences between our treatment and the original 
version of the MIT bag model. In the expression of the bag energy (1) we have 
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omitted two terms that were present in the MIT version [10,11] of the model, 



The first term was expected to represent the so-called zero-point (Casimir) 
energy. This entry was necessary to obtain a good fit in the original MIT ver- 
sion of the bag model. However, the phenomenological value of the parameter 
Zq ~ —1.9 differs substantially from its theoretical value Z ~ +0.7 [38]. Note 
that even the sign of the effect is opposite. As shown in [16], the phenomeno- 
logical value of Z can be made smaller by introducing the c.m.m. correction 
and refitting model parameters. 

The second term represents a part of the self-energy of quarks included in a 
somewhat arbitrary fashion. Such a choice reduces substantially the colour- 
electric part of the interaction energy and in the case of quarks of the same 
mass makes it vanish. This is to be contrasted to the potential model in 
which the Coulomb-like colour-electric potential plays an essential role in the 
description of the J / ip and r\ c mesons (see also Refs. [39] and [6] for the critical 
discussion on this subject). We think that the description of the same states 
in two models must not be so different. Furthermore, in the usual approach 
all the self-energy can be absorbed into the renormalization of the quark mass 
and, therefore, any use of the self-energy term in the energy expression could 
cause a double counting. So, in order to have a consistent description of the 
heavy hadrons [39] we have chosen to discard the term (22) from the bag 
model energy. 

Now let us proceed with the further modifications we want to include in our 
version of the bag model. First of all, QCD guidelines should be followed 
where possible. We will incorporate two QCD inspired modifications: scale 
dependence of the strong coupling constant a c and scale dependence of the 
effective quark mass rrif. Determination of the quark mass is a very interesting 
and complicated problem by itself. Due to the confinement, quarks are not the 
asymptotic states of QCD and, therefore, their masses cannot be measured 
directly. Moreover, the mass values depend on the chosen conventions and can 
be determined only through their influence on the properties of hadrons. The 
problem of the determination of quark masses in the context of the heavy 
quark theory is discussed in [40], and some properties of the effective quark 
mass are studied in [41]. For the recent review see the article by A.V. Manohar 




(21) 



and 




(22) 



in [42]. 
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Strictly speaking, there is no way to relate the quark mass as defined in the 
phenomenological models (such as potential model or bag model) to the pa- 
rameters of the QCD Lagrangian, or to the pole mass. Many ingredients of the 
phenomenological models are introduced by hand and can be justified only by 
the success of the model in describing the experimentally measured properties 
of the hadrons. Nevertheless, we expect these models to share some of their 
features with QCD. What can QCD tell us about the properties of the strong 
coupling constant and quark mass? From the renormalization-group analysis, 
with the rif quark flavours for which rrif « Q, in the leading logarithmic 
approximation there follows [43]: 

12% 

asiQ2) = (33 - 2n f ) M&W) ' (23) 



m 

±ln(Q 2 /A 2 ) 



MQ 2 ) = ~ ^-^r-, (24) 



where m(Q 2 ) is the mass function (running mass) in the MS scheme, A ~ 
200 MeV - the QCD constant, m - some new integration constant (analogue 
of A), and d m = 12/(33 — 2n/) - anomalous dimension of the mass. 

We are working in the soft regime where the behaviour of Eqs. (23) and (24) 
is not well-defined. So, instead of Eq. (23) we will employ the cavity-radius- 
dependent parametrization proposed in Ref. [14]: 

a - (R) = vMfhm (25) 



consistent with Eq. (23). Parameter A helps us to avoid divergences when 
R — > Rq, where Rq is the scale parameter analogous to QCD constant A in the 
momentum space. An alternative choice could be the r-dependent function ob- 
tained using the procedure adopted in [44,45] in the context of the relativized 
potential model. For the time being we prefer to use the expression (25) be- 
cause of its simplicity. We expect it to provide some average estimate of the 
scale dependence of the strength of the effective interaction inside the bag. 

For the effective quark mass we employ the parametrization 

m f (R) =m f + a s (R) ■ 5 f , (26) 



with two flavour-dependent parameters rrif and 5f. Despite rather different 
form, there is no serious contradiction between Eqs. (26) and (24). In the 
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sufficiently wide range of parameters Eq. (26) can be approximated by 



m f (R) 



771/ 



(27) 



\ ln(C/ + Rq/R) 



- d, 



with two other flavour-dependent parameters m/ and C/. Equation (27) can 
be interpreted as divergency-free extension of Eq. (24). 

Scale dependence of the quark mass proved to be important in the relativistic 
flux tube model calculations [46]. Asymptotic behaviour of the mass function 
in the framework of the Bethe-Salpeter equation coupled to the Schwinger- 
Dyson equation for the quark propagators [47] also agrees well with Eq. (24). 

It is hard to obtain equally good description of mesons and baryons in the 
quark model with the common value for the quark mass. In order to improve 
the description K. Cahill [48] proposed to use two sets of constituent quark 
masses: one set for the constituents of mesons, and another set for the con- 
stituents of baryons. Since the bag model also suffers from the flaw of this 
kind, we expect that the introduction of the scale-dependent effective mass 
would help us to improve the situation in this case as well. 

To proceed with the calculations of the hadronic properties we must relate the 
ground state energy (1) to the mass of the hadron. To this end we adopt the 
procedure proposed in [24,25] and consider the bag state \B) as a wave packet 
of the physical states \B, p) with various total momenta 



In general [49], equation of this type cannot be exact (non-relativistic harmonic 
oscillator being an exception). So, we do not expect Eq. (28) to provide the 
exact solution to the c.m.m. problem and consider this relation as a reasonable 
ansatz only. For the profile function we adopt a Gaussian parametrization 




(28) 



[26,27,32]: 




(29) 



Functions $p(s) are normalized as 




(30) 
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Parameter P that specifies the momentum distribution still needs to be deter- 
mined by making some reasonable assumption. We will use the prescription 

P 2 = lJ2^pl (31) 



where Pi = (e 2 — m 2 ) 1 is the momentum of the i-th. quark. The c.m.m. pa- 
rameter 7 will be determined in the fitting procedure. At the first sight, the 
natural choice seems to be 7 = 1, however, for the reasons that will be dis- 
cussed later we will use a more general form (31) and interpret P 2 as an 
effective momentum square. 

All averages have to be calculated with the profile $p(s). In the following we 
will need the quantities (E), (M/E), and (M 2 /E 2 ): 

(E) = J d 3 s$ 2 P (s)VM 2 + s 2 , (32) 



By using (29), Eq. (32) can be rewritten as 



(E) = J— / s 2 dsVP 2 s 2 + M 2 exp (~i 

V 7T J V 2 



(35) 



Once the energy of an individual hadron E and the effective momentum P are 
given, Eq. (35) can be solved to obtain the mass of the particle (see Ref. [20] 
for somewhat different procedure). 

As noted in [27], from Eq. (35) one can easily obtain the relation 

M 2 = (Ef-p(^)p 2 , (36) 



where 



J t 2 dtV^T^ exp (-^ 2 ) J - x. (37) 
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The limiting values of this function 0.85 and 1 correspond to ultra-relativistic 
and non-relativistic cases, respectively. Equation (36) looks much like the fa- 
miliar Einstein relation 



M 2 = E 2 — P 2 , 



(38) 



that is very popular in the various bag-model-based calculations. 



3 Results. Hadron mass spectrum 

In this section we present the calculated mass values of the ground state 
hadrons (Tables 2-4) and analyze the influence of several modifications on the 
predictions of the model. We begin with the traditional version of the MIT 
bag model [10,11]. The standard expression for the mass of the hadron in this 
model can be written as 



where the entries in the right-hand side are given by Eqs. (1), (21), and (22). 
We adopt the same model parameters B, Z , a c , m s , m c (see Table 5) as in 
the original treatment [10,11] and use the experimental mass value of the T 
meson to determine the mass of the 6-quark mj,. The up and down quarks are 
taken to be massless. 

The empirical zero-point energy term Eq used in the original version of the 
MIT bag model was later reinterpreted as representing mostly a c.m.m. cor- 
rection [16]. As the first step in modifying the model we omit this term and 
use a more elaborated procedure based on Eq. (35) to account for the c.m.m. 
In this variant of the model (denoted as Modi) the energy is given by 



which is minimized in order to determine the radius R of the spherical cavity 
in which the hadron is confined. After the minimization is performed, Eq. (35) 
must be solved numerically to obtain the mass M of the corresponding hadron. 
The free parameters of the model now are B, 7, a c , m s , m c , and m b . Instead 
of Z now we have another free parameter 7 governing the c.m.m. correction. 
To fix B, 7, and a c , the masses of the light baryons N, A, and the average 
mass of the uo-p system are employed. We use the vector meson <fi (instead 
of the baryon Q~ used in the MIT version of the model) to fix the strange 
quark mass. Our choice is motivated by an intent to have the same mass fixing 



Mmyv — E + E^jf + Eq, 



(39) 



^Modl — E + £ge lf , 



(40) 
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Table 2 

Mases (in GeV) of hadrons consisting of the u-, d-, and s-quarks in the six vari- 
ants of the bag model as described in the text. Underlined entries were used to 
define the free model parameters. x 12 was calculated without the contribution of 
the pseudoscalar mesons. 



Particle 


M cx 


Mmit 


Modi 


Mod2 


Mod3 


Mod4 


Mod5 


N 


0.939 


0.938 


0.939 


0.939 


0.939 


0.939 


0.939 


A 


1.232 


1.233 


1.232 


1.232 


1.232 


1.232 


1.232 


IT 


0.137 


0.280 


- 


- 


0.137 


0.137 


0.252 


P 


0.769 


0.783 


0.776 


0.776 


0.776 


0.776 


0.776 


U 


0.782 


0.783 


0.776 


0.776 


0.776 


0.776 


0.776 


A 


1.116 


1.104 


1.099 


1.101 


1.098 


1.116 


1.116 


E 


1.193 


1.144 


1.140 


1.143 


1.138 


1.159 


1.158 




1.318 


1.288 


1.280 


1.283 


1.277 


1.310 


1.310 


£* 


1.385 


1.382 


1.372 


1.374 


1.368 


1.388 


1.384 


"* 


1.533 


1.528 


1.513 


1.517 


1.505 


1.543 


1.536 


n- 


1.672 


1.672 


1.654 


1.660 


1.643 


1.695 


1.687 


K 


0.496 


0.496 


0.326 


0.324 


0.458 


0.437 


0.496 


K* 


0.894 


0.928 


0.896 


0.896 


0.895 


0.897 


0.895 


<t> 


1.019 


1.067 


1.019 


1.019 


1.019 


1.019 


1.019 


Xl2 




0.023 


0.021 


0.019 


0.024 


0.013 


0.012 



procedure for s-, c-, and 6-quarks. For this purpose we will employ the mass 
of the corresponding vector meson (i.e. 0, J/ip, and T). 

The experimental mass values of the hadrons were taken from the Particle 
Data Group [42] . For the isospin multiplets the averaged values were used. 

It is difficult to assess the efficacy of the different variants of the model simply 
by examining the columns of numbers presented. To assist the reader, in the 
last row of the Tables 2-4 the x N values for each of the variants are presented. 
This quantity - a root mean squared deviation from the experimental mass 
spectra - is evaluated as follows: 

i 

(41) 

where M % is the model prediction for the z-th hadron, M* x is the experimental 
value, and the sum includes N states for which sufficiently accurate values of 



X, 



1 N 

-E(m i -m: x )' 
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Table 3 

Mases (in GeV) of hadrons containing charmed quarks in the six variants of the 
bag model as described in the text. Underlined entries were used to define the free 
model parameters. 



Particle 


M ex 


-Mmit 


Modi 


Mod2 


Mod3 


Mod4 


Mod5 


A+ 

C 


2.285 


2.215 


2.259 


2.279 


2.257 


2.285 


2.285 


L. 


2.452 


2.358 


2.373 


2.393 


2.364 


2.392 


2.389 




2.469 


2.397 


2.431 


2.451 


2.429 


2.466 


2.467 


1 'c 


2.576 


2.508 


2.518 


2.538 


2.508 


2.546 


2.543 


n° 

c 


2.698 


2.654 


2.662 


2.683 


2.652 


2.696 


2.694 






3.540 


3.527 


3.552 


3.520 


3.556 


3.554 


n+ 

cc 




3.691 


3.677 


3.702 


3.671 


3.709 


3.709 


s* 

c 


2.518 


2.462 


2.464 


2.491 


2.458 


2.488 


2.482 


"* 

' 'c 


2.646 


2.604 


2.602 


2.630 


2.594 


2.637 


2.629 


c 




2.743 


2.740 


2.768 


2.729 


2.782 


2.774 


"* 

1 'cc 




3.663 


3.628 


3.668 


3.616 


3.659 


3.649 


n+* 

cc 




3.797 


3.763 


3.801 


3.751 


3.799 


3.791 


^ccc 




4.830 


4.751 


4.784 


4.738 


4.776 


4.769 




1.867 


1.726 


1.806 


1.796 


1.830 


1.833 


1.849 




2.008 


1.970 


1.994 


2.007 


1.990 


2.002 


1.998 




1.969 


1.886 


1.947 


1.936 


1.975 


1.965 


1.986 




2.112 


2.100 


2.113 


2.124 


2.112 


2.119 


2.118 


% 


2.980 


2.933 


2.999 


2.964 


3.018 


3.005 


3.020 


J/Y> 


3.097 


3.097 


3.097 


3.097 


3.097 


3.097 


3.097 






0.071 


0.042 


0.032 


0.046 


0.024 


0.027 


M ex are available. 



By comparing the results of our calculation (presented in the column de- 
noted as Modi in Tables 2-4) with the predictions of the original version 
of the MIT bag model and with the experimental values we see that for the 
hadrons containing heavy quarks (Tables 3 and 4) the agreement between pre- 
dicted and experimental values is obviously improved. For the light hadrons 
the agreement with experiment is of the quality similar to the MIT version. 
The pseudoscalar mesons (it and K) remain the source of some difficulty 
The kaon comes out about 170 MeV too low, and there is no solution for 
the pion. The small masses of these mesons result after the partial cancella- 
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Table 4 

Mases (in GeV) of the lightest hadrons containing bottom quarks in the six variants 
of the bag model as described in the text. Underlined entries were used to define free 
model parameters. The masses of B c and % were not included when determining 
X 6 . M cx of B c was taken from [2]. 



Particle 




Mmit 


Modi 


Mod2 


Mod3 


Mod4 


Mod5 




5.624 


5.548 


5.580 


5.695 


5.593 


5.624 


5.624 


s 6 


- 


5.746 


5.719 


5.839 


5.728 


5.759 


5.756 


B 


5.279 


5.148 


5.226 


5.304 


5.235 


5.252 


5.255 


B* 


5.325 


5.253 


5.283 


5.378 


5.290 


5.309 


5.306 


B s 


5.370 


5.283 


5.361 


5.435 


5.374 


5.387 


5.393 


B* s 


5.417 


5.379 


5.412 


5.504 


5.422 


5.439 


5.439 


B c 


6.32±0.06 


6.217 


6.297 


6.315 


6.307 


6.307 


6.314 


B* 




6.331 


6.335 


6.385 


6.340 


6.345 


6.345 


Vb 


9.30±0.04 


9.258 


9.438 


9.374 


9.441 


9.438 


9.442 


T 


9.460 


9.460 


9.460 


9.460 


9.460 


9.460 


9.460 






0.083 


0.033 


0.058 


0.026 


0.017 


0.018 



tion of several large terms. As a consequence, these mass values are strongly 
model-dependent and rather sensitive to the changes of the model parameters. 

Table 5 shows that the bag constant B and the mass parameters m/ have 
not changed substantially from the original MIT version. The strong coupling 
constant a c has reduced from 2.19 to 1.56. It can be noted also that for the 
light hadrons our predictions are qualitatively similar to the results obtained 
in Refs. [16,17,23]. 

Now, let us make the next step and drop out the self-energy term. The new 
version of the model (denoted as Mod2) coincides with the preceding one 
with the only exception that energy of the hadron is now given by Eq. (1) 
instead of (40). There are practically no changes in the predictions of the 
light hadron masses (only the mass of Q~ is slightly improved), while the 
model parameters undergo sizeable changes. This can be considered as some 
kind of renormalization, since the effect of the self-energy term now must be 
absorbed in the redefinition of the parameters of the model. In the heavy 
quark sector the predictions of the two versions (with and without the self- 
energy term) differ. For the hadrons with charm the overall fit is improved 
again. For the hadrons containing bottom quarks the situation is opposite. 
The agreement with experiment in the new variant is somewhat spoiled. The 
fit is still better than in the case of the original MIT version, however, it can 
hardly be considered as satisfactory. 
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Now we are in the position to examine the influence of the scale dependence of 
the effective strong coupling constant on the predictions of the bag model, by 
replacing a c with Eq. (25). In this new variant of the model we have one extra 
parameter (say, A), and therefore, in order to determine its value, some extra 
prescription is necessary. Following Ref. [29] one can put A = 1 in Eq. (25), 
and this would not be a bad choice. Another possible choice could be the 
requirement for the pion mass to vanish when m q — > [14]. Our strategy is 
somewhat different. We simply want to improve the description of the pseu- 
doscalar mesons and we can do that by adjusting the values of the parameters 
A and R that govern the behaviour of the running coupling constant (25). 
The free parameters now are B, 7, A, R , and m s , m c , m^. First, let us try to 
use the masses of four light hadrons (i.e. N, A, p-u system, and 7r) to fix the 
parameters B, 7, A, and R . Then, we employ <fi, J/ip, and T to determine 
the masses of the strange, charmed, and bottom quarks. Predictions for the 
hadron mass values generated by this version of the model are presented in 
the column denoted as Mod3 of Tables 2-4. For the light hadrons now al- 
most everything is all right. The description of the pseudoscalar mesons is 
improved considerably. The fit for the baryons is slightly worsened, but still 
remains of the quality similar to the original MIT version. The analysis of 
entries presented in Tables 3 and 4 shows that in the heavy quark sector the 
meson spectrum is improved. However, for the baryons containing charmed 
quarks the discrepancy becomes more serious. A more careful analysis shows 
that the situation may be not so bad as appears. The hadron mass differ- 
ences in the new version are described better, while the absolute position of 
the baryon spectrum is evidently positioned too low. Such regularities in the 
hadron spectrum is a welcomed feature, and we can conclude that this variant 
of the model can serve as a good starting point for a further development. 

Now we must find a way to improve the description of baryons and not to 
spoil the meson spectrum. Our proposal is to use for this purpose the scale- 
dependent effective quark mass given by the mass function (26). Instead of 
a fixed quark mass mj now we have two adjustable parameters rhf and Sf 
for each quark flavour. For fixing these parameters we employ the masses of 
corresponding vector mesons (0, J/ip, T) and the mass values of the lightest 
baryons Af containing the quark qj. The results of the fitting are given in 
the column denoted as Mod4 in Table 5. From Tables 2-4 we see that the 
agreement with experiment is improved impressively. Despite the manifest 
success in describing the hadron mass spectrum, several imperfections of the 
model still remain uncured. One drawback common to almost all variants of 
the bag model is the S-A mass difference. For the light hadrons it differs from 
the experimental value by about 30 MeV, and for the charmed hadrons the 
discrepancy of S c -A c mass splitting from the experiment grows up to 60 MeV. 
One possible solution to this problem can be the use of some chiral extension 
of the bag model [50,51,23], however, such an extension is outside the scope 
of the present investigation. 
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Table 5 

Parameters for the six variants of the bag model as described in the text. All mass 
parameters (m, in, 5) are in GeV, Ro in GeV -1 , B in GeV 4 . 



Parameter 


MIT 


Modi 


Mod2 


Mod3 


Mod4 


Mod5 


B ■ 10 4 


4.476 


4.892 


7.288 


7.597 


7.597 


7.810 


Z 


-1.836 












1 




2.480 


1.901 


1.958 


1.958 


2.004 


a c 


2.186 


1.564 


1.369 








A 








1.070 


1.070 


0.622 


Ro 








2.543 


2.543 


4.473 


m s 


0.279 


0.296 


0.347 


0.339 






m s 










0.161 


0.234 




— 


— 


— 


— 


0.156 


0.101 


m c 


1.552 


1.474 


1.614 


1.578 






m c 










1.462 


1.473 


6 C 










0.109 


0.095 


m b 


4.954 


4.696 


4.967 


4.848 






fhb 










4.786 


4.752 


s b 










0.069 


0.089 



The other problem is the high sensitivity of the mass values of the pseudo- 
scalar mesons to the changes of the parameters A and Ro that define the 
behaviour of the effective coupling constant (25). To illustrate the sensitivity 
of the model predictions upon the choice of the parameters A and R we 
present the results of an alternative calculation with somewhat different fitting 
procedure. The model parameters can be refitted to reproduce the kaon mass 
instead of the pion one. The corresponding results for the model parameters 
are presented in the column denoted as Mod5 of Table 5, and the results 
of calculations are given in the last column of Tables 2-4. We see that the 
model is rather stable in its predictions. Both versions (Mod4 and Mod5) 
provide a reasonable description of the ground state hadron spectrum. The 
most pronounced difference in the predicted hadronic mass values between 
the two versions is for the pseudoscalar mesons. It is impossible to fit the 
masses of the kaon and pion with a common set of parameters, and it teaches 
us that we cannot get all. Because the agreement with experiment is slightly 
better for the version denoted as Mod4 (with the fitted pion mass), we prefer 
to use this version of the model. However, our choice should not be taken too 
seriously. If one needs the model with the accurate kaon mass value one can 
use the version denoted as Mod5 as well. 
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There is also some concern about the masses of the heavy scalars (especially 
r]b). Light scalar mesons (77 and 77') need special treatment [10], therefore, they 
are not included in our consideration. For the attempts to solve this prob- 
lem by incorporating higher-order (~ a 2 c ) corrections see [52,53]. Meanwhile, 
it is a common practice to treat the heavy scalars on the same footing as 
the other hadrons. For the r] c meson our prediction is M(i] c ) = 3.005 GeV, 
which is about 25 MeV too high. This is an indication that the interaction 
strength for this state may be underestimated. One possible reason for this 
can be slightly too large value of the cavity radius R(r) c ) that is obtained after 
the minimization of the hadron energy. The version of the model with the 
coupling constant and quark mass fixed (Mod2) gives M(r] c ) = 2.964 GeV 
that is 15 MeV too low. By analogy we expect M(i]b) to lie somewhere in 
the region between 9.37 and 9.44 GeV (the corresponding potential model 
prediction is 9.40 GeV [45]). All these results are in some conflict with the 
recent experimental data M cx (%) = (9.30 ±0.04) GeV [54,42]. This result still 
needs additional confirmation. However, if confirmed, it could become a serious 
headache for the model builders. While the mass differences M(B*) — M(B) 
and M(B*) — M(B S ) are reproduced with good accuracy, our result for the 
775 mass value M(r] b ) = 9.44 GeV is evidently too high, signalling that scalars 
must be treated with care and in this particular case a more subtle treatment 
might be necessary. 

In order to illustrate the main features of the model, some parameters charac- 
terizing the model (for the version Mod4) are given in Tables 6-8. By inspect- 
ing the entries presented in these tables one can see how the scale-dependent 
characteristics (coupling constant, mass values of the s-, c-, 6-quarks) change 
when going from one particle to another. For example, when going from the A 
baryon to rjb meson, the strong coupling constant reduces by about 30% from 
its maximum value a max = 1.531 to the minimum value a niin = 1.046. The 
changes in the mass values are not so impressive and do not exceed 40 MeV 
for the strange, 30 MeV - for the charmed, and ~ 25 MeV - for the bottom 
quarks, respectively. 

By comparing the values of M and E one can estimate the role and size of the 
c.m.m. correction for each particle. The typical correction is ~ 400 MeV for 
the light hadrons, ~ 300 MeV for the hadrons with charm, and < 250 MeV 
for the hadrons containing bottom quarks. 

Another interesting characteristic is the function f3(x) entering Eq. (36). From 
Table 8 we see that for the bottom quarks (3(x) pa 0.99. In this case, in order 
to obtain the mass of the hadron one can use a simpler relation (38) instead of 
Eq. (35). The difference between the two values calculated using Eq. (35) and 
Eq. (38) correspondingly will not exceed 1 MeV in this case. For the hadrons 
with charmed quarks (Table 7) the difference can grow up to pa 10 MeV. This 
is not a very large difference, and the use of Eq. (38) to simplify calculations 
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Table 6 

Some characteristics of the bag model (version Mod4) for hadrons consisting of the 
u-, d-, and s-quarks. All masses and E are given in GeV, R in GeV" 1 . 



Particle 


M 


E 


R 


p(M z /P z ) 


a c (R) 


/ i i \ 

m s {R) 


N 


0.939 


1.345 


4.948 


0.925 


1.517 




A 


1.232 


1.561 


5.015 


0.944 


1.531 




7T 


0.137 


0.857 


4.417 


0.854 


1.401 






0.769 


1.155 


4.532 


0.920 


1.426 




A 


1.116 


1.559 


4.797 


0.929 


1.484 


0.393 


£ 


1.159 


1.596 


4.757 


0.931 


1.476 


0.392 




1.310 


1.784 


4.639 


0.932 


1.450 


0.388 




1.388 


1.758 


4.883 


0.944 


1.503 


0.396 


"* 


1.543 


1.950 


4.762 


0.944 


1.477 


0.392 


n- 


1.695 


2.135 


4.652 


0.945 


1.452 


0.388 


K 


0.437 


1.105 


4.189 


0.877 


1.350 


0.372 


K* 


0.896 


1.340 


4.393 


0.920 


1.395 


0.379 




1.019 


1.517 


4.271 


0.920 


1.368 


0.375 



sometimes can be justifiable in this case, too. So we see that there remains 
only the light hadron sector where the results obtained by Eqs. (35) and (38) 
may differ significantly. 

We want to end up this section with several comments. The attentive reader 
could already have noticed (look at the values of 7 in Table 5) that our effective 
momentum square (31) is about twice the value usually accepted in the various 
bag model calculations (see e.g. [16,17,22,23] etc.). It must be so, and we'll 
soon see why. 

Let us compare our method to deal with the c.m.m. with the methods used in 
Refs. [16] and [22,23]. We expect our approach to give similar results for the 
light hadrons as the others do, because the masses of the light hadrons are 
used as an input to determine the basic model parameters. In the approach 
advocated in [16] the energy of the hadron is divided into two parts: E q ~ 1/R 
associated with the quarks, and the volume energy 

4<7T 

Ey = Y BR ^ ( 42 ) 
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Table 7 

Some characteristics of the bag model (version Mod4) for hadrons with charmed 
quarks. All masses and E are given in GeV, R in GeV -1 . 



Particle 


M 


E 


R 


j3(M z /P z ) 


a c (R) 


m s (R) 


m c (R) 


A -1- 


2.285 


2.597 


4.588 


0.965 


1.439 




1.618 




2.392 


2.695 


4.572 


0.967 


1.435 




1.618 




2.466 


2.817 


4.428 


0.964 


1.403 


0.380 


1.615 


'c 


2.546 


2.885 


4.449 


0.966 


1.408 


0.381 


1.615 




2.696 


3.070 


4.336 


0.965 


1.383 


0.377 


1.612 


^cc 


3.556 


3.865 


4.205 


0.976 


1.353 




1.609 




3.709 


4.051 


4.094 


0.975 


1.328 


0.369 


1.606 


K 


2.488 


2.770 


4.668 


0.970 


1.456 




1.620 




2.637 


2.954 


4.547 


0.969 


1.429 


0.384 


1.617 


n° c * 


2.782 


3.133 


4.437 


0.967 


1.405 


0.381 


1.615 


^cc 


3.659 


3.948 


4.311 


0.978 


1.377 




1.612 




3.799 


4.119 


4.204 


0.976 


1.353 


0.372 


1.609 


CH — l~ 

"ccc 


4.776 


5.092 


3.954 


0.981 


1.296 




1.603 


D 


1.833 


2.210 


3.934 


0.953 


1.292 




1.602 


D* 


2.002 


2.327 


4.115 


0.960 


1.333 




1.607 


D s 


1.965 


2.335 


3.807 


0.951 


1.262 


0.358 


1.599 


D* 


2.119 


2.494 


3.995 


0.958 


1.306 


0.365 


1.604 




3.005 


3.397 


3.545 


0.967 


1.202 




1.593 


JH> 


3.097 


3.454 


3.689 


0.970 


1.235 




1.596 



associated with a bag. Only the part E q associated with the quarks is corrected 
for the c.m.m. 

E™=(El-P*)\ (43) 



where Pq is given by the analogue of Eq. (31) with 7 = 1. The total c.m.m.- 
corrected energy now can be written as 

E cor = E cor + Ey ^ ( 44 ) 
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Table 8 

Some characteristics of the bag model (version Mod4) for the lightest hadrons with 
bottom quarks. All masses and E are given in GeV, R in GeV" 1 . 



Particle 


M 


E 


R 




/ T~t\ 

a c {R) 


m s (R) 


/ i i \ 

m c (R) 


m b (R) 


A 


5.624 


5.777 


4.453 


0.991 


1.409 






4.883 


s fe 


5.759 


5.906 


4.478 


0.992 


1.414 






4.883 


B 


5.252 


5.425 


3.774 


0.990 


1.255 






4.873 


B* 


5.309 


5.472 


3.870 


0.990 


1.277 






4.874 


B s 


5.387 


5.591 


3.642 


0.988 


1.224 


0.352 




4.870 


B* s 


5.439 


5.632 


3.741 


0.989 


1.247 


0.356 




4.872 


B c 


6.307 


6.550 


3.297 


0.988 


1.143 




1.586 


4.865 


B* 


6.345 


6.576 


3.378 


0.989 


1.162 




1.588 


4.866 


Vb 


9.438 


9.661 


2.894 


0.992 


1.046 






4.858 


T 


9.460 


9.675 


2.950 


0.993 


1.060 






4.859 



and the minimum of this energy is assumed to be the actual hadron mass. If 
the uncorrected energy 

E = E q + E v (45) 



is minimized first (as in Refs. [22,23] and in our work), then the spurious 
energy of the center-of-mass motion is confined inside the bag, too. From the 
dimensional analysis it follows that for the massless quarks 

E = (46) 



and therefore, in order to subtract the spurious c.m.m. energy, one is forced to 
employ the relation of the type (38) in which the effective momentum square 
with 




must be used. This is an exact result (i.e. 7 = (4/3) 2 ) in a toy model with 
massless noninteracting quarks for the particular hadron chosen in the fit- 
ting procedure while determining the bag constant. In our work we treat the 
parameter 7 as a quantity to be fitted. The values obtained (see Table 5) 
favour well the Eq. (47). For the version Modi the value of this parameter 
is substantially larger because in this case we must also subtract the "spuri- 
ous" self-energy. So we see that because of its semi-phenomenological nature 
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the effective momentum square (31) cannot be associated with a pure c.m.m. 
correction, but it may also contain some other ~ 1/R corrections. 

We have established the link between the treatment of Ref. [16] and ours. Now 
it is evident that the version Modi of our treatment is similar to the model 
used in [16] and, therefore, in the case of the light hadrons, both models should 
give qualitatively similar results. The prescription advocated in Ref. [16] seems 
to be somewhat more physical. However, this approach is hardly compatible 
with the Eq. (35) we have used to employ the c.m.m. correction. 

The approach adopted by authors of Refs. [22,23] is more similar to ours. 
They minimized the "c.m.m.-uncorrected" energy with the self-energy term in- 
cluded, and used Eq. (38) to incorporate the c.m.m. correction. However, their 
"uncorrected" energy contains the so-called zero-point energy term, which can 
be interpreted as representing mostly a c.m.m. correction [16]. Their zero-point 
constant Z ps —0.8 is much smaller in comparison with the original MIT value 
Z —1.84. So, the expression of the energy adopted in Refs. [22,23] can be 
considered to be partially c.m.m.-corrected. In other words, in their approach 
the c.m.m. correction is incorporated in two steps. First, the term Z /R is 
subtracted from the energy and this partially corrected energy is minimized. 
Then, Eq. (38) is applied, in which the usual expression for the momentum 
square of the quarks confined in the bag, Pq = J2 n iPh * s used. Eventually, 
one obtains similar results as in Ref. [16]. At present we have no simple an- 
swer to the question what quantity, the energy or the mass (c.m.m.-corrected 
energy), must be minimized (see the discussion on this subject in Ref. [19]). 
As we have seen, in practice this is somewhat a matter of taste, and seemingly 
rather different methods may give quite similar results. 



4 Electroweak properties 



The bag model also sets a framework to calculate other static properties of 
the hadrons. In this section we present our results for some electroweak prop- 
erties: magnetic moments, charge radii, and axial-vector coupling constant. 
To obtain some feeling for the sensitivity of computed quantities to the model 
assumptions, we compare our predictions with the results of the original MIT 
model and with the experiment. Magnetic moments of the hadrons in the 
static spherical cavity approximation can be represented in the form [9] (see 
also [32]): 

^ = EK(h|^|h), (48) 
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Table 9 

Composition of baryon magnetic moments. The label I is used to collectively repre- 
sent up and down quarks. 



Particle h 


< 

'11 


P 




N 


-y i 


£+ 


±(8^ + //) 


S° 


h(2fi l +n s ) 


£~ 


\u s - V) 


A 




H° 








S° A 




Q- 




A++ 


111 1 



where is the charge of the i-th quark and parameters // are given by [10]: 

i _ 4:EjR h + 2mjR h - 3 R h 
~ 2(e i R h -l)e i R h + m i R h 6 ' 1 ' 



In the last expression £j represents the energy of the i-th quark and i?h stands 
for the bag radius of the hadron under consideration. Magnetic transition 
moments are defined by 

/4Uh' = (hVk* |h> , (50) 



where R^ = R^ is assumed. 

Matrix elements (h'| a % z qi |h) can be calculated with SU(6) wave functions as 
described in Ref. [55], and for the cases we are interested in are displayed in 
Table 9. 

The square charge radius and axial- vector coupling constant can be calculated 
from the expressions 

Rh 

\ " f 2 



I 



^ / r 2 dr Pf( r ) + Ql(r) , (51) 



o 
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9°A 



dr 



(52) 



Analytic expressions for Eqs. (51) and (52) can be found in Ref. [10]. In the 
case of massless quarks the value of g° A does not depend on the radius and 
is always equal to 1.088. 

Before comparing the static quantities with the corresponding experimental 
values they must be corrected for the center-of-mass motion. For this pur- 
pose we adopt the prescription proposed in Refs. [24,25] (see also [32]). Their 
formulae for the corrected values of /i, (r 2 ), and Qa are [24]: 



1 + (M/E) + {M 2 /E 2 } 
3 

h ~ 1 + 2{M 2 /E 2 ) 
3 



• o 1- (M/E) M P 

H o ~Vh 



M h 



9gh 

4 p2 



9a 



1 + 2 (M/E) 



9a- 



(53) 
(54) 
(55) 



In the equations above, the values of (M/E) and (M 2 /E 2 ) are defined by 
Eqs. (33) and (34), respectively, P 2 is given by Eq. (31), M P is the mass of the 
proton, Q h and M h stand for the charge and the mass of the corresponding 
hadron. 

Our predictions for the magnetic moments are presented in Table 10. In order 
to simplify the calculation of the transition moment /i(S° — > A), the same wave 
function (that of the £ baryon) was used for both states. The experimental 
values were taken from the Particle Data Group [42]. 

Predictions for other electroweak properties are given in Table 11. The exper- 
imental values for rp, r^_, r n , and tk were taken from Refs. [56], [57], [58], 
and [59], respectively. 



Predicted values in the variants Mod2-Mod5 of the model are of similar quality 
(differences between the calculated values do not exceed 5%), therefore, we list 
only the values obtained using variants Modi and Mod4. Agreement with the 
experimental values is good in both variants of the model. Predictions for the 
electroweak parameters calculated in the variant Mod4 are comparable to the 
results obtained in Ref. [22]. 

Predicted values of the electroweak parameters practically are insensitive to 
the corrections associated with the scale dependence of the effective coupling 
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Table 10 

Magnetic moments of baryons (in nuclear magnetons). Uncorrected values are en- 
closed in parentheses. 



Particles 




MIT 


Modi 


Mod4 


P 


2.79 


1.90 


(2.08) 


2.89 


(1.88) 


2.61 


N 


-1.91 


-1.26 


(-1.39) 


-1.85 


(-1.25) 


-1.66 


£+ 


2.46 


1.83 


(2.01) 


2.67 


(1.75) 


2.34 






0.58 


(0.63) 


0.81 


(0.54) 


0.70 


XT 


-1.16 


-0.67 


(-0.75) 


-1.05 


(-0.66) 


-0.94 




-1.25 


-1.05 


(-1.13) 


-1.44 


(-0.95) 


-1.22 




-0.65 


-0.43 


(-0.45) 


-0.64 


(-0.36) 


-0.54 


A 


-0.61 


-0.48 


(-0.51) 


-0.66 


(-0.43) 


-0.55 


S°^A 


-1.61 


-1.08 


(-1.19) 


-1.54 


(-1.05) 


-1.35 


fi- 


-2.02 


-1.54 


(-1.56) 


-1.91 


(-1.26) 


-1.57 


A++ 


3.7-^7.5 


4.16 


(4.36) 


5.40 


(3.81) 


4.76 



Table 11 

Some electroweak parameters of the hadrons. All charge radii are given in fm. Un- 
corrected values are enclosed in parentheses. 



Parameter 


M ex 


MIT 


Modi 


Mod4 


r P 


0.88 


0.73 


(0.79) 


0.80 


(0.71) 0.70 


r 


0.77 


0.68 


(0.76) 


0.77 


(0.66) 0.66 




0.66 


0.49 






(0.63) 0.66 


r K 


0.58 


0.47 


(0.65) 


0.70 


(0.58) 0.60 


9A 


1.26 


1.09 


(1.09) 


1.34 


(1.09) 1.33 



constant and the quark mass. In contrast, the c.m.m. corrections for mag- 
netic moments and for axial- vector coupling constant improve the predictions 
significantly. The corrections for the charge radii in all cases are of minor im- 
portance. Neutron charge radii in our version of the model remain zeroes. This 
drawback of the model is a direct consequence of the isospin symmetry. 

Our feeling is that one must not take the good agreement of the corrected 
values of the electroweak parameters with the experiment too seriously. For 
example, the version Modi, in which the c.m.m. correction seems to be over- 
estimated, better agrees with the experiment than the version Mod4. In fact, 
because we have employed the c.m.m. correction in a somewhat phenomeno- 
logical fashion, we cannot disentangle the pure center-of-mass motion and 
other possible effects. The treatment of c.m.m. in our work is far from being 
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perfect. Moreover, other effects, such as recoil corrections and the pion cloud 
contribution, may also be important. Nevertheless, despite of all this criti- 
cism, the model seems to provide reasonable predictions for the electroweak 
properties of the hadrons. 



5 Summary 

One of the objectives of this paper was to examine the influence of the correc- 
tions associated with the center-of-mass motion, the scale dependence of the 
running coupling constant, and the scale dependence of effective quark mass 
on the mass spectrum and on other static properties of the hadrons, calculated 
in the framework of the bag model. Special attention is paid to the hadrons 
containing heavy (charmed and bottom) quarks. All quarks are treated on 
equal footing. The heavy quarks rattle inside the bag cavity in the manner 
of the light ones with the maximum of their distribution being closer to the 
center of the bag than the same of the light quarks. 

Incorporating the corrections consecutively we were able to investigate the ef- 
fect of these corrections upon the predictions of the model. We have found that 
the proper treatment of the center-of-mass motion is essential to obtain the 
reasonable description of the heavy hadrons. The running coupling constant 
and scale-dependent effective quark mass proved to be useful ingredients of the 
modified bag model. These corrections helped us to obtain the good agreement 
of the calculated masses of the heavy hadrons with the experimental values. 
There is strong evidence that these two modifications of the model should be 
applied simultaneously. 

The bag model with all these corrections included can be treated as rather 
simple and controllable framework for the unified description of the light and 
heavy hadrons. Maybe the worst discrepancy of the model is the tt-K mass 
difference. Another systematic discrepancy inherited from the original MIT 
version of the model is the Xh~A h mass splitting. Finally, the description of 
the r/ h states also seems to be somewhat problematic. Despite the several draw- 
backs mentioned above, the model accounts reasonably well for the masses of 
almost all hadrons under investigation. The accuracy achieved in the descrip- 
tion of the hadron spectrum suggests that for the further improvement an 
explicit breaking of the isospin symmetry may be necessary. 
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